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Abstract 

The possibility of baryogenesis at the electroweak phase transition is considered 
within the context of a minimal supersymmetric standard model with spon- 
taneous R-parity violation. Provided that at least one of the sneutrino fields 
acquires a large enough vacuum expectation value, a sufficient baryon asym- 
metry can be created. Compared to R-parity conserving models the choice of 
soft supersymmetry breaking parameters is less restricted. The observed baryon 
asymmetry, ^f- ~ 10~ 10 , can be explained by this scenario and the produced 
baryon-to-entropy ratio may easily be as high as — ~ 10~ 9 . 
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One of the major open questions in theoretical particle physics is to explain the 
observed baryon asymmetry of the universe. The well-known Sakharov conditions 
for baryogenesis, namely T, B, C and CP -violation |T[, can be fulfilled in the early 
universe during a phase transition for suitable models of particle physics. In particular, 
it has been shown that the Standard Model (SM) includes all the requirements for 
electroweak baryogenesis 0. However, the electroweak phase transition is too weakly 
first order to preserve the generated baryon asymmetry in the SM ||. Furthermore, 
the CP-violation due to the phase of the Cobayashi-Maskawa -matrix is too small 
for a sufficient generation of baryons . The realization of electroweak baryogenesis 
therefore requires physics beyond the SM. 

The Minimal Supersymmetric Standard Model(MSSM) has appeared to be a 
promising candidate for the explanation of the origin of the observed baryon asymme- 
try in the universe. It has been shown that with an appropriate choice of parameters 
a large enough baryon asymmetry, — ~ 1CT 10 , could have been created at the elec- 
troweak phase transition [[| || . The necessary choice of parameters involves, however, 
quite stringent bounds. Furthermore the choice of |sin0 M |, where <f)^ is the phase 
of one of the soft supersymmetry breaking parameters //, has been constrained even 
more in a more recent analysis of the bubble wall profile 0. It is therefore useful and 
interesting to expand the idea of electroweak baryogenesis to less constrained super- 
symmetric models. In this paper we shall consider the creation of baryon asymmetry 
in the context of a model where R-parity is spontaneously violated ||. 

The Minimal Supersymmetric Standard Model (MSSM) assumes the conservation 
of a discrete symmetry, R-parity || [U|. R-parity is related to the spin (S), total 
lepton (L) and baryon (B) number of a particle by R p = (— i^ 3B + L + 2S so that all the 
standard model particles have an even R-parity whereas their supersymmetric partners 
have an odd R-parity. There is, however, no a priori reason to expect that R-parity is 
conserved. Therefore it can be a spontaneously or even explicitly broken symmetry. A 
model that spontaneously violates R-parity proposed in ref. JTT|] has a superpotential 
of the form (suppressing all SU(2) and generation indices) 

h u QH u u c + h d H d Qd c + hJH d e c + (h H u H d - /i 2 )$ + hjH u 9 c + h$S9 c + h.c. , (1) 



which conserves both total lepton number and R-parity (hat denotes a superfield). 
The new superfields ($, Of, Si) are SU(2) <g> U(l) -singlets and carry conserved lepton 
numbers (0, —1, —1) respectively. 

In this paper we shall focus on the h e lH d e c term. Using standard methods |9j , one 
can show that the following terms are present in the Lagrangean: 



£r = H[h yi ji>* Li R + h ek jV* Rk L]ej } 



(2) 



where 

ii I/, ; (3) 




and R, L are the chirality projection operators. 

The Higgsino current associated with charged higgsinos can be written as || 

J% = H-fH. (4) 

The Higgsino current (^) is associated with a triangle diagram where Higgsinos interact 
with electron-like leptons, similar to the diagram in ref. 0. The CP- violating source 
in the diffusion equation is in this case created by interactions arising from the term 

C R = h eij Hv*Lej, (5) 

where Vi =< Ur& >. It is worth noting that a non- vanishing source requires only 
that one of the sneutrino fields acquires a vev. This is due to a complex vev so 
that two degrees of freedom are present. In a spontaneously R-parity breaking model 
the global B-L symmetry is spontaneously broken by assuming non-zero vev's for 
Vr =< vrz >, vg =< S3 > |TTj]. In general we may also choose a non-zero value for 



vl =< vis >■ However, spontaneous R-parity breaking generates a majoron that are 
produced in a stellar environment and to suppress the stellar energy loss one must 
require ^ m w ^. 10~ 7 JTT[] . This is easily achieved for vr = Oil TeV) provided that 
vl < 0(100 MeV). For generality we consider here the case where all left- and right 
handed sneutrinos have a non-zero vev. In numerical estimations we make the usual 
assumption that only vrs is large. 

If one chooses interactions similar to ref. |J where both left- and right-handed elec- 
trons are included the amount of CP- violation is proportional to vrPl -term. However, 
if we consider only interactions like (|5|), the total CP- violation has a factor u R , which, 
assuming that \vr\ ^> \vl\, can give a much larger contribution. 

The scenario considered here sets some limitation to the hierarchy of the elec- 
troweak phase transition, namely that the higgsino and sneutrino field(s) acquire their 
respective vev's simultaneously. Otherwise sphaleron transitions wash out the previ- 
ously created baryon asymmetry. 

The CP-violating source may now be computed using CTP-formalism as described 



in ref. [12]. The interactions of the Higgsinos lead to a contribution to the self-energy 



of the form 

^cp{ x i V) = 9cp( x > y)RG° i> ' < (x, y)L + gcp(x, y)LG° i> ' < (x, y)R, (6) 



where 



3c p = KiM x ) h ekjV* k (y). 



(7) 



Similarly for the other component High temperature corrections change the 



dispersion relation of charginos and neutralinos [ 13] . The spectral function of Higgsinos 
may, in the approximation r^Cm^, be written as [12] 



p s (k,k°)=i^+mfj) Uk +ie+iY fl f-u J l{k)) - ((k -ie-iT 6 ) 2 -oj%(k) 



-i 



-i 



where oJ 2 ^{k) = k 2 + m^(T). The effective squared Higgsino plasma mass m 2 g(T) may 
be approximated by m 2 ^(T) « \fj,\ 2 . Similarly, should be replaced by the electron 
plasma mass -^g 2 T 2 for p ej (k, k°). 
The source term 

f-T 



S S = -fd 3 r 3 f tr[Yl> P (X,x 3 )G%(x 3 ,X)-G%(X,x 3 )^ P (x 3 ,X) 



-oo 

1> 



+ G%{X, x 3 )E^ P (x 3 , X) - £< P (X, x 3 )G%(x 3 , X)] 
now contains the following function 

fg = 9cp( X i x s) - 9cp(x3,X) 

= «Im [(^^^(^(X)^^) - ^(x 3 K(X))], 

which after the Higgs insertion expansion |j| becomes 

/, = flm [(^^)(^(X)^^(X) - ^(X)d^-(X))]. 



(9) 



(10) 



(11) 



The possibility of complex Yukawa coupling constants, h e , deserves further discussion. 
The h e matrix clearly has 18 degrees of freedom i.e. 9 complex phases are present. 
From (B|) we note that 7 of the complex phases may be set zero by redefining the 
sneutrino, higgsino and electron-like lepton fields. Generally we then have two complex 
phases left. For simplicity, we shall from now on assume that all Yukawa coupling 
constants are real i.e. h ei j G M, i,j = 1,2,3. However, the sneutrino vev's may be 
complex, so that we may write 



i6i{X) 



Vi {X) = Ai{X)e 

where Ai(X),9i(X) e 1R. Function f g can now be written as 

f g = i(h eij h ekj )\sm(9i - 9 k )[AA' k - A k A[] - cos(0; - 9 k )AiA k (9i + 9 k )' 



(12) 



(13) 



The first term is of the form presented in || with an extra factor of sin(#j — 0k)- The 
second term is of a new form, which exists due to the non-vanishing complex phases 
of the sneutrino vev's. 

We can now write the expression for the CP-violating source, 



Sfj = 2i(h eij h ekj ) [sm(0* - k )[AiA' k - A k A[] - cos(0, - k )AiA k (0i + 9 k )']l^, (14) 



where XA is slightly modified from the result in [HJ, 

k 2 (u} e jUjj — k 2 ) 



r * = L°° dk 2%^* [('-a.W)^^,r, 

+ (l - 2Re (nH))F(u e .,r ej ,u 67 V%) + 2(Jm (n#) + Im (n e 

£(^ r #,^,r e .)l ( 15 ) 



where n k = [exp (^ ( ,/T + ir#, J + 1] and 



F(a,b,c,d) = ^[(a + c) 2 + (6 + cf) 2 ] _1 sin[2 arctan 

+ \\( a ~ c ) 2 + ( b + <i) 2 ] _1 sin[2 arctan ^-j-^] 



G(a,b,c,d) = — -[(a + c) 2 + (b + d) 2 } 1 cos [2 arctan 7— ~r] 

- \[{a - c) 2 + (b + d) 2 ]- 1 cos[2 arctan ^— (16) 

The damping rate of Higgsinos and electron-like leptons is dominated by weak inter- 
actions so that we may take « r e ~ 0.05 T 0. 

We shall now review the basic ingredients of electroweak baryogenesis as presented 
in ref. ||. In the paper the baryon to entropy ratio was shown to be 

s v*s 

where g(ki) is a numerical coefficient depending on the degrees of freedom, D the 
effective diffusion rate, r„, s = 6k«^T (k = 1) Jl4| the weak sphaleron rate and v w the 



velocity of the bubble wall. The entropy density s is given by 

2n 2 g* s T 3 



s 



45 



where g* s is the effective number of the relativistic degrees of freedom. A was shown 
to be a weighted integral over the CP-violating source 7(1*) = v w f(ki)d u J°(u): 

1 r°° 

A= ~mT du^u)e~ x + u , (19) 
D\+ Jo 



where A + = (v w + yv^ + 4TD) / (2D) and the wall was defined to begin at u — 0, 
where u denotes the co-moving coordinate u = z + v w t and the wall is assumed to 
move in the direction of the positive z-axis. (f(ki) is a coefficient depending on the 
number of degrees of freedom present in thermal path and related to the definition of 
the effective source || |T]|.) In the notation of the current paper and after performing 



a partial integration (assuming the CP-source vanishes at u = and u — > oo) the 
coefficient A is related to the shape of the bubble wall through 



/•CO 

A oc 7i = / duSfr{u)e~ 
Jo 



2i(h eij h 



ekj, 



du 



sin(^-0 fc )[A^-A fc ^] 



cos(0i - fc )>Mfc(0i + O k )']e- x+u )l?l. 



(20) 



In comparison, in the paper by Carena et al. |6| the coefficient A is related to the 
shape of the bubble wall through 



A oc h 



oo P) 

du—iH.H^ - H 2 H[)e- x +\ 
ou 



r°° d 
/ du-—(H(u 
Jo ou 



du 



(21) 



where if = (if^ + H 2 )^, tan/? = Hi/H 2 and iij's are the real parts of the neutral 
components of the Higgs doublets. In the paper by Riotto the CP- violating source 
associated with the charginos was considered and in that case I 2 can be written as 



A oc I 2 = Im (n) [ / duH{uf^-e- x ^\ [3M 2 g 2 2 lf + M 1 g 2 1 l§] 



(22) 



where Mi are the soft supersymmetry breaking parameters, X^' B are corresponding 
integrals for the winos and the bino similar to (|15|). 

(|20|) can be written as (writing out the sums explicitly) 



,e ~ x+u du 



j i 

/■oo 

+ E M*i / t sin (^ - e k)(AA' k - A k A[) 



cos(^ - ^)AA fc (^ + e k )']e- x + u du 



-T e 3 



= ^^(Sx + ^ + ^r 



(23) 



In the wall shape was assumed to take a simple sinusoidal form where the field 
H{u) can be given by 



H(u) 



1 — cos 



un 



[6{u)-6{u-L w )]+v6{u-L w ) 



(24) 



and the angle [3{u) by 
R( \ ^ 



1 f u7r 

1 — cos I — 



[6{u) - 6{u - L w )\ + Ap9(u - L v 



(25) 



where A/3 is given by Aj3 = (3{Tq) — arctan(mi(T )/m2(To)), calculated at the tempera- 
ture where the curvature of the one-loop effective potential vanishes at the origin. (It is 
thus the angle between the flat direction and vacuum direction.) We shall assume that 
both the modulus and the phase of the sneutrino vev's have a similar sinusoidal form 
In this approximation, one can write an analytical expression for the integral in 



vr 2 (2 X% Li + 20 XI Lj vr 2 + 3 (ll + 5 ^) ^ ^ 
2 K 4e A + L » (A^ L^ + 14At 7r 2 + 49A 2 h L 2 u 7r 4 + 36 7r 6 ) ^ 

= h{L w )A^v\ (26) 

where L w is the width of the bubble wall and v 2 =< iJi > 2 + < H 2 > 2 . 



Similarly we can now estimate the first term in fl23|) . 



^i = -E^ 2 %/i(^), (27) 

i 

where 9f i =< ^ >. Here it is assumed that the phases of sneutrino fields acquire their 
vev's sinusoidally from < 0, >= to < 6>j >= Of.. £> 2 ,3 cannot be solved analytically. 
However, before numerical analysis, some special cases are worth a further study. Let 
us define new phase variables, i[)i k (X) = 0i(X) — k (X) and (f)(X) ik = 0i(X) + 9k(X). 
Now, if ipik{X) = ipik, Bx t 3 clearly vanish and we are left with an expression similar to 
( PS] ) with an additional sinf/^ factor, 

/+ = + Ai E heijh ekj fi(L w ) sin ^ ik {v] + v 2 k )A(3 ik X^ , (28) 

j k<i 

where Af3 ik is the angle between the flat and the vacuum direction at the origin. If, 
on the other hand, ipi k {X) = i.e. all sneutrino fields acquire the vev of their phases 
in uniform, we are left with 



If = -4*E [Y, h lij v i e fh{L w ) + 2'£hei j h ekj v i v k e f f 1 (L w )]rJ, (29) 

j i k<i 

where Of =< 0i > (all phases have the same vev). Clearly since B 2 is the only positive 
term and sin^ < 1, 1^ gives the maximum value of I\. Similarly If gives the 
minimum value. However, the overall sign of I± is not significant so that we are only 
concerned with the absolute value of I\. Since If and If have opposite signs, one of 
them give the maximum value of |ii | depending on the choice of parameters. Generally 
If is greater in magnitude than If due to the A/3 ik factor in If. 

To compare the magnitude of the CP-violating source considered in the present 
paper with previous results, numerical estimates are in order. In choosing parameters 
we shall follow closely ref. |J. The bubble wall width is chosen to be L = 25 /T (our 
results are quite insensitive to the exact choice of the bubble wall width). 
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Figure 1: Values of 1^ at different phase transition temperatures. 



First we shall consider (H). We choose M 2 ps ps Mi ps T ps fOO GeV and 
11 ~ Xj to give us an estimate of the effect of the CP-violating source. With these 
values sin^ is roughly of the order 0.1 and M 2 can numerically be shown 

to have a value of about 10. Putting all these values into (J^), we get an estimate 
I 2 ~ 0.35f 2 A/3. Estimating v = 246 GeV and A/5 ps 0.015 we finally arrive at an 
estimate I 2 ~ 300. Recent calculations using a more realistic bubble wall shape have 
lead to an additional suppression factor of ~ 0.3 [ff|. This effect is ignored here since 
the analytical approximation is used in both cases. 

Similarly we can now estimate (|23"D . Choosing L w = 25 /T we first find that 
fi(L w ) ps 0.25 for T ~ 100 GeV. Values of ig are shown in fig. [l]. The values 
are very indifferent to the lepton type considered so that we may very accurately take 



= Ifj. As from fig. 1 can be seen, if we choose \i ps T 



100 GeV, Xfj is 



about —2. From the figure it is clear that there exists a resonance at value /x ps T. 
However, this resonance is not as pronounced as in ref. ||. 

Numerical examinations show that the integral in B 2 is approximately linear in ipik 
and we may estimate it with ~ 0.18^ (no significant T dependence). Furthermore the 
integral in B% is approximately constant and we may take it to be ~ 0.05. Substituting 
these values into ft 



we then have 



k<i 



h « -i8 [ - °- 25 E h % v l e fi + °- 18 E h eiMM + v l) &Mik - 0-05 £ Vi v k <f> ik 

k<i 

(30) 



Here we have several unknown parameters. To estimate the numerical value of (|30|), 
let us first assume that v 3 ^> vi,v 2 and 6^ = 9f 2 ~ i.e. only the tau-sneutrino (we 



may have as well chosen any other generation) has a significantly large vev. In this 
approximation 



h « -i%v\O h [ - 0.25 h% + 0.18(A/3 3 i £ ^sAu + A/3 32 £ /^A?;)] • (31) 

3 3 3 

If we now choose v 3 = 1000 GeV, estimate h ei j ~ 10 _1 and take a conservative estimate 
for Afiij ~ 0.1, we get an estimate l\ ~ —50000 9f s , which is significantly greater 
than the previous estimate provided that 9f 3 > 0.01. The values of the different 
contributions are in this case If ~ 10000 sm9f 3 and I{" ~ —60000 9f. A . Even if 
h ei j ~ O(10 -2 ) and 9f 3 ~ 0(1), I\ is still large enough to explain the origin of the 
baryon asymmetry of the universe (even when accounting for the extra suppression 
factor due to a more realistic bubble wall profile). 

We can now easily estimate the possible baryon to entropy ratio resulting from the 
spontaneously R-parity violating model. In || it was estimated that with the same 
choice of parameters used in this paper, — ~ 10~ n . Since I1/I2 can easily be ~ 10 2 
the baryon asymmetry created within the context of the scenario considered here can 
be as high as ~ 10~ 9 . So the required value, y ~4x 10 -11 , is accessible with a non- 
restrictive choice of parameters. 

In the present paper we have estimated the contribution to CP-violation in the 
bubble wall at the electroweak phase transition arising from the terms present in min- 
imally supersymmetric spontaneously R-parity violating models. It has been shown 
that the amount of CP-violation due to the higgsino-lepton -interaction can quite 
easily be larger than the contribution due to the higgsino-gaugino -interaction by a 
factor of 10 2 . This contribution is due to the possible existence of a large (complex) 
sneutrino vev and a Yukawa coupling constant of order O(10 _1 ). It is notable that 
this result is not dependent on the soft supersymmetry breaking parameters and 
furthermore allows freedom in the choice of the value of /^-parameter which may also 
be chosen real. Also the phase transition temperature dependence is quite weak so it 
can be concluded that this mechanism allows for the creation of baryon asymmetry at 
the electroweak phase transition with quite a flexible choice of parameters. However, 
it should also be noted that there are several sources of uncertainty present due to the 
approximate nature of the chosen bubble wall profile and how higher corrections to 
the CP-violating source behave. Furthermore, the weak sphaleron rate is under dis- 
TBj and changes on that may significantly change the results on electroweak 



cussion 
baryogenesis 
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